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Abstract. - We determine the maximal work extractable via a cyclic Hamiltonian process from 
a positive-temperature (T > 0) microcanonical state of a N 3> 1 spin bath. The work is much 
smaller than the total energy of the bath, but can be still much larger than the energy of a single 
bath spin, e.g. it can scale as 0(V N In N). Qualitatively same results are obtained for those cases, 
where the canonical state is unstable (e.g., due to a negative specific heat) and the microcanonical 
state is the only description of equilibrium. For a system coupled to a microcanonical bath the 
concept of free energy does not generally apply, since such a system — starting from the canonical 
equilibrium density matrix pr at the bath temperature T — can enhance the work extracted from 
the microcanonical bath without changing its state pr- This is impossible for any system coupled 
to a canonical thermal bath due to the relation between the maximal work and free energy. But 
the concept of free energy still applies for a sufficiently large T. Here we find a compact expression 
for the microcanonical free-energy and show that in contrast to the canonical case it contains a 
linear entropy instead of the von Neumann entropy. 



How much work can be extracted from a state of a phys- 
ical system via cyclic processes? This question governs our 
understanding of energy conversion and storage, and hence 
is central for thermodynamics [I}JS]. The basic answer, 
known as the Thomson's formulation of the second law, is 
that an equilibrium state cannot yield work. This formula- 
tion is an axiom in thermodynamics, but its first-principle 
derivations were given in literature for a canonical (Gibb- 
sian) equilibrium state 2 . The main consequence of the 
Thomson's formulation is that only non-equilibrium states 
can be sources of work. The maximal work extractable 
from such states via a cyclic process was studied both for 
macroscopic [T] and finite systems [1H5]. 

One instance of the maximal work is especially well- 
known, because it provides the physical meaning of free 
energy [T] . Consider a quantum system with Hamiltonian 
H and initial density matrix p. This system is in contact 
with a canonical thermal bath at temperature T. External 
fields act cyclically on the system + bath. Assuming no 
system-bath coupling both initially and finally, the maxi- 
mal work extracted by the fields reads [TJ[H] 



F[p] - F[p eq ], F[p]=tv(pH) + TtT(p]np), (1) 



where F[p] is the free energy and p oq = e 



- H / T /tr[e~ H / T ] 



is the canonical equilibrium state of the system, which 
is its final state after work-extraction pQ. The maximal 
work is determined by the deviation of p from its canonical 
equilibrium value p cq as quantified by the free energy (JT|). 

One notes however that all above results refer to a spe- 
cific notion of equilibrium, viz. the canonical state. An- 
other concept of equilibrium is given by the microcanoni- 
cal state, which describes an isolated system equilibrated 
due to its internal mechanism [TU] , or an open system cou- 
pled weakly to its environment (so weak that no energy 
is exchanged) [5J[TT]. This is a more fundamental notion 
of equilibrium: (i) under certain conditions the canonical 
state can be derived from it for a weakly coupled sub- 
system pp. (ii) In contrast to the canonical state, whose 
preparation refers to an external thermal bath, the micro- 
canonical state can be applied to a closed few-body sys- 
tem provided that it satisfies certain chaoticity features 
[T2l[13] - (Hi) Since local stability conditions of the canoni- 
cal state are more demanding — a fact closely related to the 
no work-extraction feature [I] — there are situations, where 
the equilibrium can be described by the microcanonical 
state only, since the canonical state for them is unstable 
[14] . For such systems, frequently realized via long-range 
interactions, the entropy is a non-concave function of en- 
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ergy, and hence the notorious macroscopic equivalence be- 
tween canonical and microcanonical state is broken j!4j . 
Even if this macroscopic equivalence holds, it is by no 
means obvious that in the argument around (JT|) one can 
substitute the canonical state of the thermal bath by the 
microcanonical state [with the same temperature] , because 
in general the work (pp) is not a macroscopic quantity, i.e. 
it does not scale with the number of bath particles. This is 
however widely done in literature, e.g., when introducing 
the free energy as in (pp) one basically never specifies the 
equilibrium state of the bath; see, e.g. [UE]. 

We revisit the maximal work-extraction problem for a 
thermal bath in a quantum microcanonical state. It was 
noted already some work can be extracted via a cyclic 
Hamiltonian process from a few-particle microcanonical 
system [3] . Recent papers studied to which extent the ex- 
traction of work from one-particle classical microcanonical 
system can be carried out by physically realistic Hamilto- 
nians [15(116] . Our purposes here are different: 

• We focus on finding the maximal amount of work 
extractable from a macroscopic microcanonical state of 
N S> 1 particle thermal bath. 

• We also determine the work extracted via a sys- 
tem coupled to a microcanonical thermal bath, and check 
whether the reasoning (JT} generalizes at least qualitatively, 
i.e. whether the concept of free energy applies to the mi- 
crocanonical situation. 

The subject of work extraction via a system coupled to 
a thermal bath is an active research topic. Refs. [16-20 
discuss various set-ups for this problem: quantum, clas- 
sical, with or without state-dependent feedback etc. Re- 
call that {1} is at the core of relations between statistical 
thermodynamics and information theory [7]. The term 
tr([p — p C q]H) in the right-hand-side of (pp) is the energy 
extracted from the system, while the remaining (entropic) 
part comes from the bath. If tr([p — p cq ]H) is negligible 
(e.g., because H contains only few almost degenerate en- 
ergy levels), the work is extracted from the bath due to 
the difference between the initial entropy and its canonical 
equilibrium value. This relation between the entropy and 
work is the essential part of information driven engines 
(e.g., Szilard's engine) [7]. In contrast, various forms of 
fuel operate due to the initial non-equilibrium energy, i.e. 
the term tr([p — p C q]H) in (pp). 

Microcanonical thermal bath. The microcanonical 
state is characterized by two parameters: energy E and 
width a pp. The corresponding density matrix is diago- 
nal in the energy representation, all energies within the 
interval [E, E — a] have equal probabilities, all other en- 
ergies have probability zero. For a N ^> 1-particle system 
the number d(E, a) of energy levels within the interval 
[E, E — a] defines the microcanonical entropy [1] : 

S(E) = \nd(E,a) = 0(N), (2) 

where the choice of a should not influence the leading 
O(N) behavior of S(E). Eq. © is the von Neumann 



entropy for the microcanonical density matrix ^ ; see also 
[2Tj . For clarity we want to work with a specific model of 
a macroscopic microcanonical system (bath). This is the 
basic model of the field: N ^ 1 uncoupled two-level spins; 
each spin has energies and 5 > [8] . Some of our results 
extend to more general bath models, as seen below. 

The bath Hamiltonian reads (diag[. . .] means diagonal 
matrix in the energy representation) 

H = diag[0, S, 26..., 5N], (3) 

where each element Sk is repeated dk times, 

N\ 

dk 3 k\{N-k)V (4) 
Hence every energy shell 8k is dfc-degenerate. Denote 

e fc = (l,...,l), fc = (0,...,0), D k =Y] k d m , (5) 

^ — ^m— 

where (0^) is the vector of k l's (O's). 

For the present model of bath the microcanonical state 
is easy to define: all energies SM have equal probability 
all other energies have zero probability. Thus we put 
a — > 0, the minimal thermodynamically consistent width 
for this model. Note that the degenerace of the energy 
levels is convenient, since it allows to set a — > 0. It is 
however not essential: an effective degeneracy will be any- 
how regained for a small but finite a > 0, since the energy 
levels of a macroscopic system are located very densely [I] . 
The bath initial state reads in representation ([3]) 

= -i- diag[0D M _ 1 ,e dM) Djv - £ , M ]. (6) 

For N ^> 1 this microcanonical state does have desired 
features expected from thermodynamics, e.g. macroscopic 
equivalence with the canonical state, equilibration of a 
small subsystem, third law 18 . The density matrix of a 
single bath spin is Gibbsian ex diag[l, e~ s / T ] with [8] 

e - s/T = m/(l-m), m — M/N. (7) 

The same T is recovered as microcanonical temperature 

m 

l/T = dS(E)/dE, E = SM, S(E)=lnd M , (8) 

where S(E) is the microcanonical entropy This equiv- 
alence can be shown via formula (jTTJ) that is proven below. 

Note that although the spins are uncoupled, the micro- 
canonical state does not reduce to the tensor product of 
the separate spin states (otherwise it would amount to 
the canonical state). It contains inter-spin correlations. 
Ultimately, this is the reason why, as seen below, a micro- 
canonical bath can yield work in a cyclic process. 

We restrict ourselves with M/N < 1/2, i.e. positive 
temperatures. The case with M/N > 1/2 is definitely less 
interesting, because now each spin of the bath is in a state 
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with a negative temperature. Such states are trivially ac- 
tive, i.e. they yield work in a cyclic process. 

Work extraction. At some initial time t — the 
bath Hamiltonian H{t) becomes time-dependent due to 
interaction with sources of work. Consider a cyclic process 



H(0) = H(t) = H, 



(9) 



where r is the final time. The work extracted in this 
thermally isolated cyclic Hamiltonian process is 

W = tr(i?[fii - fi f ]) = 5M - tr(HQ{), Q £ = U£liU\ (10) 

where S7f is the final state of the bath, and U = 
q~ e -(' l / h ) Jo j- means chronologization. Conversely, 

for a given unitary U one can construct a class of Hamil- 
tonians that generate U and satisfies ((9]) [5]. 

Condition is necessary for the system to be an au- 
tonomous carrier of energy that should deliver work to 
another system (e.g., to a work-source) via an interaction 
which switches on and off at well-defined times. Hence 
this is a cyclic Hamiltonian process. 

We now maximize the work W — or minimize the final 
energy tr(7Jf2f) — over all cyclic Hamiltonians, i.e. over 
unitary operators U. Note from (flOl) that 

tr(JTfif) = V 2 1 EaCabibl^b), Cab = \(b\U\a)\ 2 , (11) 

z — 'a, b=l 

2 N o N 

C ab >0, V C ab = V C ab = l, (12) 

^ — J a— 1 — ^o—l 

where {E a }a=i and {|a)}„ =1 are, respecitively, the eigen- 
values and eigenvectors of H [see ((3])], and the elements 

Three conditions in (fT2l) 



are defined in 

mean that the matrix C a b is double-stochastic [55]. Con- 
versely, every such matrix can be represented as C a b = 
\(b\U\a)\ 2 for some unitary U Every double-stochastic 
matrix equals to a convex sum of permutation matrices 
nH (Birkhoff's theorem [22]): C = ^A^H, £ q A q = 
1) A a > 0, where each matrix nl a l acting on a column- 
vector x amounts to permuting (in a certain way) the 
elements of x. Eq. (fTTj) shows that tr(iJfif) is a linear 
function of the matrix C — {C a b}- Hence its minimum 
over the unitary operators U, that is its minimum over 
double-stochastic matrices Cab, is reached for C ab equal 
to some permutation matrix n. It is clear from ([lip that 
n, when acting on the vector (6|Oi|6) permutes its ele- 
ments such that all its non-zero [equal to each other] ele- 
ments concentrate at lowest energies {E a }a=i 0- For the 
final state we have (a|Of |a) = ^ b Il a 6(6|f2i|6). Hence the 
lowest-energy final state compatible with i7; reads 



Of = -7— diag[e dju ,0 Dj , 



(13) 



Once n is found we can employ the standard procedure 
of constructing the corresponding unitary operator U and 
the cyclic Hamiltonian 5 . Note that n does depend on 



the energy of the initial state fi;: n applied on a micro- 
canonical state with a different energy will not lead to the 
maximal work-extraction. This does not differ from (say) 
the ordinary Carnot cycle, whose implementation also de- 
mands knowing the initial state of the working body. 
The maximal work W max reads from ([TBI [TU|) : W ma x — 



[El=o ekd k + (M~£+l)(d 



5M- 



summation by parts, 



M ' 



Et'odk)]. After 



1 



1 ^M-e 



Di 



(14) 



where is defined in ([S]), and where integer I = l(M) is 
found from 



dk > d M > 2^ „ d k . 



Eq. (HH [in]) hold for any microcanonical state (0 [5]) ; the 
specific form (£[]) is not necessary. 

We shall now calculate W max f° r t wo limits: T — > oo 
and a finite N, and then N — > oo and a finite T. 

Doubly maximized work. We set the number of 
spins N to a large, but a finite number, and maxi- 
mize W max (T') over all positive temperatures of the N- 
spin bath. The maximum is reached for T = oo (or 
M = N/2 as shows) and provides an upper bound for 

the work extractable from the positive temperature bath. 

— -t 

We now calculate W mllx (oo). Consider the sum^ fe 2 =0 dk = 

N_ 

J2m=edN._ m in (fT5|) . The dominant summation region is 
m ~ I. We shall see below that I <C N. Hence for ((3]) 
we use the Gaussian approximation dN__ m = 2 2 e N ' 2 
[23] . change the sum to integral, and find £ from 1 = 



JV 



In 



JV 



41n(JV/4) 



N 



o 



lnln(jV/4) 
ln(Af/4) 



J^dxe'Wm-. £ = 

The second term under square root is negligible if In N is 
large. Eq. (HU) then implies (for N > 1): 



xxe W2J 



SI: 



-VAT In AT. (16) 



This is a reachable upper bound for the work extractable 
from the A^-spin microcanonical bath. 

Finite temperatures and thermodynamic limit. 

We employ ([4]) and assume the standard thermodynamic 
limit: m = M/N < 1/2 (and hence T > 0) in is a fixed 
finite number for M,N — > 00. We note from and 
that for any fixed finite numbers m, £ and A^ — > 00, 



i -^=e-^[i+o(h; 

CtjVm A 



(17) 



Since the sums in (fT4l H"5j) are dominated by their largest 
terms, using (|17l) (with m < 1 and integer £) amounts to 
calculating these sums via geometrical progression, e.g., 
£Et~o +1 dk = ElZi-l e~ 5k/T - We get for £ and W max 



W max (T) = 5[£-l + v e (1-v)- 

-S/T p = \ H l - v ) 

In v 



(18) 
(19) 
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where \x] is the ceiling (upper) integer part of x, e.g., 
[0.99] = 1, [-0.99] = 0. According to Q3JJ), £ grows to 
infinity with T: £ = 1 for e - & ' T < |, £ = 2 for > 



-5/T > 



3 for 0.68232 > e^/ T > 



i+\/5 



Wmax(7 1 ) is a continuous function of T, but 



etc. 
W 



has 



jumps at the temperatures, where £ changes, and one en- 
ergy shell in the final density matrix (fT3"|) is completely 
filled; see (fT5"]) and Fig. Q] Hence 

W max (T) = 0(5) for N > 1 and T = 0(<5). (20) 

The situation is symmetric with respect to different spins 
of the bath. Hence after the work-extraction the initial 
energy of each bath spin changes negligibly = 0{-^). The 
final state of each spin is diagonal in the energy represen- 
tation and thus after work-extraction it has a well-defined 
temperature that differs from the initial temperature by 
O ( jt ) . Recall that the Gibbsian state as such cannot yield 
work in a cyclic process © [2J. Hence the work is ex- 
tracted due to inter-spin correlations present initially in 
the microcanonical state. 

Eq. (fT5|) shows that W max (T) increases faster than T: 



W max (T) = T 







Hi) 


+ 1 







for T > 6, 



(21) 



where we used 



ln(l- 



ln 1 



\n(l-v) 
In v 



Eq. (|2"Tj) is practically 



good already for T > 1.5 5. Now W max can be much larger 
than the energy of a single bath spin. In the limit T 3> 5 
this energy is equal to 6/2; see ([7]). 

Microcanonical states not equivalent to the 
canonical one. Eq. ([20| does not depend on the con- 
crete form (pEJ) of (ife. What is needed for (|2"0"|) is that the 
sum ^2ff =0 dk is dominated by its last term cLm- Then 
(fl5]) implies £ = 0(1), and CH} leads to ([20]). Hence ([20]) 
generalizes the Thomson's formulation of the second law 
to the microcanonical situation. In particular, ([20]) holds 
for those dk, where the macroscopic equivalence between 
microcanonical and canonical states is violated. As an ex- 
ample consider (J3][6]) with d,M = e. N ^ M ' N > . This spectrum 
satisfies all above conditions and leads to (|20[) . Now the 
entropy lnc?M is a convex function of energy 5M. Hence 

T i-i 



the specific heat C 



l (SM) i 



calculated from 



is neg- 



ative, and the macroscopic equivalence between canoni- 
cal and microcanonical states is clearly violated, because 
C > is an automatic consequence of the canonical state 
PP. Such convex-entropy spectra are realized in macro- 
scopic long-range interacting systems |14) . 

Another example of convex entropy and canonical- 
microcanonical non-equivalence, where still (|20l) holds, 
is the first-order microcanonical phase transition [Tlll4j. 



where in the vicinity of some critical energy E, 

S(E) 



S(£) 



has 



a jump: ^-j^Ib 



E 



tJt. This de- 

i i 



scribes coexistence of two phases with different temper- 
atures. Since a more stable phase should have a larger 



entropy, we get <T\ [T] . The above non-equivalence is 
seen here, because at a canonical first-order phase transi- 
tion different phases have the same temperature pp . Even 
though two phases at different temperatures do co-exist, 
the extracted work has the same order of magnitude (|20l) 
as for a homogeneous-temperature microcanonical state. 

System coupled to the bath. An important instance 
of the maximal work problem is the amount of work ex- 
tractable from a thermal bath in the presence of a smaller 
system coupled to it; sec (pp). How much work can be 
extracted from a combined state of a two level system 
with energies and e > and the microcanonical thermal 
bath? Answering this question will alow us to understand 
to which extent the concept of the free energy applies to 
the microcanonical situation. Before starting the analysis 
we should like to stress again that so far the statistical 
physics literature does not distinguish between canonical 
and microcanonical situations when introducing and ap- 
plying the free energy concept; see e.g. pp. 

Let the initial density matrix of the two-level system be 
Pi] its eigenvalus are ttq > The spectrum of the overall 
initial state i?i = pi <g> f2i reads [see © O] 

1 

d M 



Spec[i?i] = — [7roe dM ,7rie dM ,02D 



-2d» 



7T >7ri. (22) 



Both initially and finally the two-level system and bath do 
not interact. Hence the overall Hamiltonian % reads 



n 



H(t) =H s ®l + l(g)H 
diag[0, S, 28 . . . , 8N, e, e + 5, . 



(23) 

,e + SN], (24) 

where H is given by ([3]) , and H$ is the two-level Hamilto- 
nian with energies and e. We recall that each symbol kS 
(or e + kS) in ([24]) is repeated dk times. Once we consider 
unitary work-extraction processes, the final state of the 
overall system will have the same eigenvalues (1221) . Re- 
calling our discussion between (|TT|) and (THfl) it should be 
clear that the minimal final energy for the overall system 
is achieved for the unitary operator that forces Rf to have 
the same eigenvectors as H and permutes the eigenvalues 
(|2"2l) such that the largest eigenvalue is matched with the 
smallest energy, next to the largest eigenvalue with the 
next to the smallest energy and so on. Note that Spec[i?i] 
is already ordered in a non-increasing way. It remains to 
order 124|) in a non-decreasing way and write the lowest 
final average energy as scalar product of two vectors 

tr(URi) = Spec[i?i] • [0, 5, . . . , a5, e, (25) 
{a + 1)5, e + 6,..., 5N, e + 6(N - a),. . . ,e + 6N], 

where a = [_f J , and \_x\ is the floor (lower) integer part of 
x, e.g., L0.99J = 0, L-0.99J = -1. 

Consider the work extracted from the overall system 
that is maximized over all unitary dynamic operators. 
Since the system and bath do not interact both initially 
and finally, this work separates into two parts coming, re- 
spectively, from the system and bath [see (ITU1 IT4"]) ]: 

tr (n [Ri -Rf]) = W max + W SUI +tr(H s [ Pi - Pi ]), (26) 
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Fig. 1: Thick curve: VK max given by (|14p as a function of tem- 
perature T for 5=1. Normal curves show Want given by ([27] 
as a function of T for 8 = 1, 7Ti = 0.3 and (from top to bottom) 
a= LfJ =0, 1, 2. 



where tr(iJg[pi — pf]) is the energy change of the two- level 
system. W max + W SUI is the work coming from the bath. 
Here W max is given by (|14|) (the maximal work extracted 
from the bath alone) and we defined the surplus work W S m 
(work extracted from the bath, but due to the system) . 

Obviously, W sur + tr(_ff s [pi — Pf]) > 0, since tr(H[i?; — 
Rf] ) results from optimizing over a larger set of parameters 
than W max . Note that tr(_ffg[pi — pf]) appears also in the 
right-hand-side of ([T]), and there is some analogy between 
Wsm and the entropy difference T(tr[— p eq lnp eq + plnp]) 
in which is the work extracted from the canonical 
bath. There the work W max extracted from the canonical 
equilibrium bath alone (without the system) is zero. 

Scalar product (|25|) is calculated straightforwardly; for 
clarity we focus on the thermodynamic limit regime f|17[> : 



%iF 2 + (7ro-7n)Fi], (27) 
e[n 1 P 2 + {n Q -Ti 1 )P 1 l (28) 

where W max is given by (fT8|) and we defined for k — 1, 2: 



ur 

tr(tf s p f ) 



1 - i; 



1 



1 - v l + v° 



+asign(^i fe - £ 2 k) 

£ lk +a 



v'™ (1 + v 1+ot ) 



Pk = 



V 



1 - V 
In 



+ sign(^i fe - t 2 k) 

( k(l-v 
\ 1+v" 



k- 



1-v 



(29) 



1-v 



In v 



i2 k = 



In ( fc d-") 
In v 



(30) 



Recall that sign(0) = 0, a = |_4 J , and that [ccj and 
are defined after (j2"5"j) and (|T5)) . respectively. 

The final state pf of the two-level system is diagonal in 
its energy representation. The eigenvalues of pf are read- 
off from (|28|) . The excited state of pi is less populated 
than the ground-state; otherwise it can still provide work 
via a cyclic process. In this specific sense the two-level 
system partially equilibrates; recall that p; is arbitrary. 

Fig. Q] displays W max (T) and W sur (T) for a represen- 
tative range of parameters. It is seen that for a = 0, 



Wgur > and both W max and W sur monotonically increase 
with T. For a > the positivity of W SUI is recovered only 
for a sufficiently high T provided that ir ^ m ; see Fig. [TJ 
For 7r = 7Ti, we always get W SUI < 0. 

Consider now e < S [i.e. a = 0in (|25p] and assume 
V = e _5 / T < 1/3 for simplicity. Eqs. (I27H5D)) produce 



W sur = 
tr(FsPf) = 



vS(2ir Q -l)(l-v)- 2 , 
e[»7T + (l-3«)7ri](l 



«)" 



(31) 
(32) 



Eq. (f3~Tj) shows that, in addition to W / sur +tr(iJs[A~Pt]) ^ 
0, the work extracted from the bath is enhanced, M^ sur > 0, 
for any state of the two-level system besides the completely 
mixed one, where 7To = 7Ti = \. 

The energy difference tr(Hs[pi — Pf\) can be positive or 
negative. Hence parameters can be tuned such that it is 
zero, e.g., from (|3"2"]l and for e = Tin 2 < Tin 3 < 5 we 
get that initially canonical equilibrium two-level system, 
7To = 1 — 7Ti = (1 + e _e / T ) _1 = |, enhances the work ex- 
tracted from the bath without changing its marginal state: 
tr(ifgpf) = 0. Hence for enhancing the work extracted 
from the microcanonical bath one needs that the system 
is ordered: its state should not be completely mixed, while 
the maximal enhancing is achieved for a pure state. But 
the state of the system need not change. 

A non-equilibrium system coupled to canonical equilib- 
rium bath can enhance the work extracted from the bath 
only at the cost of changing (towards equilibrium) its ini- 
tially non-equilibrium state; see ([lj. The free energy mea- 
sures this change. For the microcanonical bath, the work 
can be enhanced already by an equilibrium two-state sys- 
tem without changing its marginal state. We conclude 
that the concept of the free energy does not generally ap- 
ply to a system coupled to a microcanonical bath. 

But this concept applies in the high-temperature limit. 
For T > 5, e we get from ([27H30]) and from flST): 



W SUI = (l-27ri)rin2, 
tv(Hspi) = e/2. 



(33) 
(34) 



Eq. (|34|) means that the final state of the two-level sys- 
tem is completely mixed, which for the present high- 
temperature case coincides with the canonical equilibrium 
state. Eq. (|33p predicts work-enhancement only if initially 
the two-level system was out of equilibrium [recall that 
\ > tti]. Hence for T>iJ,ewe recover the logics of the 
canonical-bath situation, but not its letter, because for a 
canonical bath H^ sur reduces to the difference between two 
von Neumann entropies that are logarithmic functions of 
the initial eigenvalues ttq and tti ; cf . our remark after (|2"rj|) . 

We shall show elsewhere that a for a /x-level system in a 
state (density matrix) p with eigenvalues ordered as ttq > 
7Ti > ... > 7T^_ i, we can define the linear entropy as 



£[p]=y M 1 7rJ(fc + 1) mfifc + 1) -kink}. 

*k= 1 



(35) 



Generalizing (|33|) . the surplus work W SVLT extracted from 
a high-temperature microcanonical bath in contact with 
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this system is then W sur = T(ln/x - C[p]) = T(£[|] - 

C[p\), where i is the maximally mixed state of the p- 
level system. For the canonical situation this expression 
involves the von Neumann entropy — tr[plnp] instead of 
C[p\. Note that ln^t > C[p] > 0: the upper (lower) limit 
is reached for the maximally mixed (pure) p. 

Hence for a system in initial state p and Hamilto- 
nian H$ coupled to the high-temperature microcanoni- 
cal bath one can define the microcanonical free energy 
F[p] = tr(Hsp) — T£[p], whose difference F[p] — ^[pcq] 
(after adding to W max extracted from the bath alone) de- 
fines the maximal work extracted from the system+bath. 

Summary. We reformulated the Thomson's formula- 
tion of the second law for a JV > 1 particle equilibrium 
bath in a microcanonical state: if the bath temperature 
T is finite, the maximal work extractable from the bath 
via a cyclic Hamiltonian process is > <5, where S is the 
energy of a single bath particle. The maximal work tends 
to Sy/N In N if N is large but fixed and T oc. The 
reformulation applies equally well to both ordinary mi- 
crocanonical states, which are macroscopically equivalent 
to canonical states, and convex-entropy microcanonical 
states for which no canonical state can be defined, e.g., 
because of a negative specific heat [14]. The existence 
of such states demonstrates that a viewpoint on a micro- 
canonical state as emerging from measuring the energy of 
the canonical state is not generally valid. Thermodynam- 
ics of such systems can have peculiarities [24] , but we saw 
that they satisfy the same generalized Thomson's formu- 
lation much in the same way as ordinary microcanonical 
states. The work extraction is possible, since the micro- 
canonic state of the bath is not Gibbsian, though each its 
constituent can be in a Gibbsian state. 

It is widely known that only a non-equilibrium sys- 
tem can lead — at expense of changing its state towards 
equilibration — to work extraction from a canonical bath 
PQ. This work is given by the free energy difference (fT]). 
In contrast, a canonical equilibrium system (having the 
same temperature as the bath) can enhance the work ex- 
tracted from the microcanonical bath without changing its 
marginal state. Hence the concept of free energy, in the 
sense of the maximal work, does not generally apply to the 
microcanonical situation. The application of the concept 
is recovered for T ^> 5, but the canonical expression of the 
free energy is not restored, instead it should be formulated 
via the linear entropy (|35[) . 

* * * 
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